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Dephasing and relaxation of the nuclear spins coupled to the nitrogen-vacancy (NV) center during optical 
initialization and readout is an important issue for various applications of this hybrid quantum register. Here 
we present both an analytical description and a numerical simulation for this process, which agree reasonably 
with the experimental measurements. For the NV center under cyclic optical transition, our analytical formula 
not only provide a clear physics picture, but also allows controlling the nuclear spin dissipation by tuning an 
external magnetic field. For more general optical pumping, our analytical formula reveals significant contribu¬ 
tion to the nuclear spin dissipation due to electron random hopping into/out of the m = 0 (or m = ±1) subspace. 
This contribution is not suppressed even under saturated optical pumping and/or vanishing magnetic field, thus 
providing a possible solution to the puzzling observation of nuclear spin dephasing in zero perpendicular mag¬ 
netic field [M. V. G. Dutt et al., Science 316, 1312 (2007)]. It also implies that enhancing the degree of spin 
polarization of the nitrogen-vacancy center can reduce the effect of optical induced nuclear spin dissipation. 
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I. INTRODUCTION 

Diamond nitrogen-vacancy (NV) center— is a leading plat¬ 
form for various quantum technologies such as quantum com¬ 
munication, quantum computation, and nanoscale sensing^rl. 
The electronic spin of the NV center and a few surround¬ 
ing nuclear spins form a hybrid quantum registes^—. Impor¬ 
tant advantages of this solid-state quantum register include the 
long electron and nuclear spin coherence timsi^, the capabil¬ 
ity of high-fidelity initialization, coherent manipulation, and 
projective readout of the electronic/nuclear spins^ii^ and even 
the entire quantum registe r*h'"^d^ by optical and microwave 
(or radio frequency) illumination. However, during the optical 
illumination for initialization and readout^iSiiii^^, the dissi¬ 
pative spontaneous emission and non-radiative decay of the 
NV electron generates substantial noise on the nuclear spin 
qubits through the hyperfine interaction (HFI), which may 
significantly degrade the control precesion. This motivates 
widespread interest in using the NV center electron to engi¬ 
neer the nuclear spin dissipation, including pure dephasing 
and relaxationi^iii. 

In the past few years, the optically induced nuclear spin dis¬ 
sipation has been investigated in many work a^^'^^i*^~ —. Gen¬ 
erally, the nuclear spin dissipation originates from the ran¬ 
dom fluctuation of the NV electron under optical illumination, 
which falls into two categories; one involving the flip of the 
NV electron spin and the other does not. The former is usually 
strongly suppressed by the large energy splitting of the NV 
electron unless the NV electron is tuned to the ground state or 
excited state level anticrossin g^^'^* . The latter is energetically 
more favorable and dominates the nuclear spin dissipation in 
many situations, as confirmed by a series of experiment a^'^^’ . 
The theoretical investigation of this latter mechanism has been 
carried out in the framework of a phenomenological spin- 
fluctuator modeU^. This work gives an intuitive understanding 


for the optically induced nuclear spin dissipation; the gener¬ 
ation of a rapidly fluctuating effective magnetic field on the 
nuclear spins by the optically induced random hopping of the 
electron between different states. When the hopping is suffi¬ 
ciently fast and hence the noise correlation time is sufficiently 
short, the nuclear spin dissipation could be suppressedi^ in a 
way similar to the motional narrowing effect in NMR spec¬ 
troscopy in liquids. This effect has been successfully used to 
significantly increase the nuclear spin coherence timei^. 

Despite these remarkable success, this spin-fluctuator 
model still suffers from two drawbacks. First, its analytical 
form is qualitative, while obtaining quantitative results require 
numerical simulations. This not only complicates the calcula¬ 
tion, but also smears the underlying physics picture. Second, 
the various parameters in this model are phenomenological, 
i.e., they are not directly related to the physical parameters of 
the NV center, but instead are obtained from fitting the ex¬ 
perimental data. This precludes a straightforward guidance 
on controlling the nuclear spin dissipation by tuning various 
experimental parameters. 

To bridge this gap between experimental observation and 
theoretical understanding, we present a microscopic and an¬ 
alytical theory on the nuclear spin dephasing and relaxation 
by an optically illuminated NV center at room temperature. 
In addition to performing numerical simulation of the cou¬ 
pled NV-nuclear spin evolution, we further derive analytically 
a closed Lindblad master equation for the nuclear spin by adi- 
abatically eliminating the fast electron spin dynamics in the 
Born-Markovian approximation. We begin with the simplest 
case in which a single cyclic transitions (e.g., between the 
ground and excited m - Q states) of the NV center is optically 
driven. Our analytical expressions for the nuclear spin de¬ 
phasing and relaxation provide a quantitative description and 
a physically transparent interpretation that substantiates the 
previous analytical (but qualitative) and numerical resultsi^. 
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They also demonstrate the possibility to control the nuclear 
spin dissipation by tuning the magnetic field^. Next we con¬ 
sider general optical illumination of the NV center incorporat¬ 
ing finite inter-system crossing between m - 0 and m = +1 
subspaces. Our numerical results agree well with the experi¬ 
mental measurementsiii. Our analytical results shows that the 
random hopping between the m = 0(orm = +l) triplet states 
and the metastable singlet of the NV center could significantly 
contribute to nuclear spin dissipation. This contribution is not 
suppressed under saturated optical pumping and is nearly in¬ 
dependent of the magnetic field. This provides a possible so¬ 
lution to the puzzling observation of nuclear spin dephasing 
in zero magnetic field^. An analytical formula for the nuclear 
spin dissipation in terms of the HFl tensors also allows us to 
measure the HFl tensor for the excited electron state, which 
is usually smeared by the short electron spontaneous emission 
lifetime. 


II. TWO-LEVEL ELUCTUATOR MODEL: ANALYTICAL 
RESULTS 

A. Model 


To begin with, we present a microscopic theory for the de¬ 
coherence of an arbitrary nuclear spin I (e.g., '^N, or '^N) 

by the electron of the NV center undergoing optically induced 
cyclic transition |g) |e), e.g., |g) = |0) and \e) — \Ey) in 

the widely used setup for single-shot readoutii22. In the rotat¬ 
ing frame, the electron dynamics is governed by the Liouville 
superoperator XeCO = -i[He,{-)] + ZaJaS^iLaX-), where 


D.R 

He - ^&e,e + 


is the electron Hamiltonian, it/j s A is the detuning 

of the optical pumping, and is the rate of the crth dissipa¬ 
tion process La in the Lindblad form D{La){-) = La(-)Ll - 
(■))/2. Here we include the spontaneous emission 
L = &g e from \e) to |g) with rate y\ » 1/(12 ns) and the pure 
dephasing L - (Je,e of the excited state \e) with rate y,^, which 
has a strong temperature dependence, from a few tens of MHz 
at low temperature up to 10^ MHz at room temperature^^. 
Including the electron-nuclear HFl (S^Ag-i-SeAf.)! s Fland 
the nuclear spin Zeeman term y^B ■ I (y^r = -10.705 kHz/mT 
is the nuclear gyromagnetic ratio) under a magnetic field 
B, the electron-nuclear coupled system obeys 

p^LeP-i[(t + yN^)-lp] ( 1 ) 

in the rotating frame of the pumping laser. 

There are two contributions to the nuclear spin dissipa¬ 
tion. One involves the flip of the electron spin and hence is 
strongly suppressed by the large electron-nuclear energy mis¬ 
match away from the NV center ground state and excited state 
anticrossings. The other does not flip the electron spin and 
hence is energetically favorable in most situations. In our 
analytical derivation, we neglect the former contribution by 
dropping the off-diagonal electron spin flip terms in F and 


only keep the diagonal part: F « ^g,g<^g + ^e,e<i^e, where 
(Dg - (gISgIg) ■ and oje - {e\Se\e) ■ Ag. The second term 
of Eq. ([T]i describes the precession of the nuclear spin with 
angular frequency yA^B -H a>g and yA^B -H cjg, respectively, con¬ 
ditioned on the electron state being |g) and \e). When a>g + (Og, 
the optically induced random hopping of the electron between 
Ig) and \e) gives rise to random fluctuation of the nuclear spin 
precession frequency and hence nuclear spin dissipation: the 
fluctuation of the precession frequency orientation (magni¬ 
tude) leads to nuclear spin relaxation (pure dephasing)i^. Be¬ 
low we derive analytical a closed equation of motion of the 
nuclear spin to describe these effects. 


B. Lindblad master equation for nuclear spin 

The time scale for the optically pumped two-level NV cen¬ 
ter to reach its steady state is ~ tnv = 1/(2^ -i- yi) < 12 ns, 
where/? = 27r(QR/2)^i5*-^'>'‘'^'>''"^'^^^(A) is the optical pumping rate 
from Ig) to |e) and 5^^'(x) = (y/7r)/(x^ + y^) is the broadened 
^-function. When tnv is much shorter than the time scale of 
the nuclear spin dissipation, we can regard the NV center as 
always in its steady state P as determined by J^gP - 0, e.g., 
the steady state population on \e) and |g) are Pg - RI{2R + yi) 
and Pg - I - Pg, respectively. Then we treat the dissipative 
NV center as a Markovian bath^l and use Born-Markovian 
approximation to derive a Lindblad master equation for the 
reduced density matrix of the nuclear spin p{t) = Tig p(t) (see 
appendix A for details): 

p = -i{(b ■\,p]+ 2F^2)[/z]p + T+D{I+\p + F_D[/_]p, (2) 

where cb = yAiB -H PgCOg + PgiOg is the average precession fre¬ 
quency that defines the nuclear spin quantization axis Cz = 
tt)/|tj|. The last three terms describe the nuclear spin dissipa¬ 
tion in the tilted cartesian coordinate 


Cx = Cj: sin (,0 - Cj, cos (, 0 , (3a) 

Cf = cos ip cos 9e.x + sin p cos 6e.y - sin flC;, (3b) 

Cz = (bl\<b\ - sin 0cos + sin 0 sin + cosflC;, (3c) 

where 0 {p) is the polar (azimuth) anlge of cb in the conven¬ 
tional coordinate (e^:, e^, e^) with e, along the N-V symmetry 
axis. The nuclear spin dissipation include pure dephasing [the 
second term of Eq. (|2])] due to the fluctuation of Pz and relax¬ 
ation [the last two terms of Eq. (|2]), with I± = Ix ± ily] due to 
the fluctuation of P± = Px ± iPy- Typically the nuclear spin 
level splitting |u)| <K yi, y^, so we obtain 

Tg 2 

E^o = “ (^g)z\ , (4a) 

T-2 

Ty^T-^^\{cOg-cOg)A\ (4b) 

where = Ox^x + Oy^y is the component perpendicular to 
the nuclear spin quantization axis, T = 1 //? H- l/(yi H- /?) is the 
duration of one electron hopping cycle (excitation time 1 //? 
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and de-excitation time 1 /(yi + R)), and 


^ ^ + 7ry2^«r.+r.)/2)(A) ^ ^ 

R + yi 2R + yi 2R + yi 

( 5 ) 

is the uncertainty of the electron dwell time in the excited state 
in each hopping cycle. Here the last step of Eq. Q holds at 
room temperature, where y^ ~ 10^ MHz is much larger than 
typical yi,R, and A. Equation (|4| shows that nuclear spin 
dissipation vanishes when cOg = (Og, simply because in this 
case the nuclear spin precession frequency is not randomized 
by the optically induced electron hopping. 


C. Physical picture 

Equations (|2l)-(|5]l not only provide an quantitative and an¬ 
alytical description for the dissipative nuclear spin dynamics 
due to an optically pumped NV center, but also have a phys¬ 
ically transparent interpretation that substantiates the previ¬ 
ous analytical (but qualitative) and numerical resultsi^. Eor 
example, the pure dephasing rate in Eq. (l4al l is directly con¬ 
nected to the nuclear spin phase diffusion process by the op¬ 
tically induced random hopping of the electron between the 
ground state |g) and the excited state \e)^. To clearly see this, 
lets consider the phase accumulation of the nuclear spin dur¬ 
ing an interval [0, f]. Suppose that during this interval, the 
electron undergoes N hopping cycles, and that during the kth 
cycle, the electron stays in |g) for an interval t*:, so the to¬ 
tal dwell times in |g) and |e) are t = and t - t, re¬ 

spectively, and the nuclear spin accumulates a phase factor 
g-/(aj+yA,B)zT-i(a,+yA,B)z(»-T) ^ ^ number of hopping 

cycle ss f/r » 1 , i.e., t is the sum of many random vari¬ 
ables {tic}, so t obeys Gaussian distribution centered at Pgt 
with a standard deviation VfVrg, where Tg is the rms fluctua¬ 
tion of each Averaging the phase factor over this Gaussian 
distribution gives where E^ coincides with Eq. (l4al l 

as long as Tg is given in Eq. (|5]), e.g., at room temperature, for 
weak pumping R yi, the uncertainty Tg =» V 2/71 of the 
dwell time in \e) is dominated by the uncertainty in the spon¬ 
taneous emission; while for saturated pumping, Tg =» 1 /(V 2 /?) 
is strongly suppressed by the rapid optically induced transition 
between \e) and |g). The relaxation rate r± in Eq. ( l4b] l can be 
understood in a similar way. 

analytical results Eqs. © provide a microscopic basis 
for the previous modeli^ and experimental observation a^i*^'^^ , 
e.g., it clearly shows the initial increase of the dissipation rates 
r^, r± cc R under weak pumping R ^y\ and the motional nar¬ 
rowing r^,r± oc 1 //? under saturated pumping R :s> y\. The 
former arises from the increase of T with decreasing R under 
weak pumping, while the latter comes from both the decrease 
of Tg ~ 1 /R and T ~ IjR under saturated pumping. Our ana¬ 
lytical formula also demonstrate the possibility^^ to control Tg, 
and r± by using the magnetic field to tune the nuclear quan¬ 
tization axis Cz oc u>, e.g., if we tune Cz to be perpendicular 
(parallel) to o)g — a>g, then we can eliminate nuclear spin pure 


dephasing (relaxation). Interestingly, the sum rule 

2 

+ = ( 6 ) 

suggests that reducing E^ (r±) inevitably increases r± (E^j,) 
and it is impossible to suppress E^ and r± simultaneously, un¬ 
less the NV states are tuned such that oig - a>g. 

D. Connection to experimental observations 

Equation (|2|l describes the dissipative evolution of the nu¬ 
clear spin in the tilted coordinate (ex, Cy, Cz) with ez oc u>. 
Erom Eq. I©, we obtain the Bloch equations 

(h) 

= (7a) 

11 

d,{h} = m - ^)<4), (7b) 

j 2 

for the average nuclear spin {1(f)) = Trlp(f), where Ti = 
l/(r+ -i-r_) and T 2 - l/(ri^ + (r+ H-r_)/2). Then the sum rule 
in Eq. (|6]l implies I/T 2 + l/(27’i) oc \{i>g - i.e., tuning 

the magnetic held can prolong T 1 time (T 2 time) at the cost of 
reducing T 2 time {Ti time). 

The above Bloch equations have simple solutions = 

{/z(0))e“7A and{/+(f)) = However, nuclear 

spin initialization and measurement are usually performed in 
the conventional coordinate (e^,, e^, e^) with along the N-V 
axis, so Ti and T 2 will be mixed in the observed signals. Eor 
example, Dutt et al? initialize a strongly coupled nuclear 
spin-1/2 (hereafter referred to as '^Cb, according to the nota¬ 
tion of Gali^i) into the eigenstate (| t) + I i))/ V2 of 4, let 
it evolve freely for an interval t, and then measure through 
a 7r/2 pulse followed by a Huorescence readout of 

via the NV center. According to Eq. O), the measured signal 
{hit)) - Y^a=X, yz(e.t ■ ^aXIait)) consists of a non-oscillatory 
term gjjj 2 ^^2 that decays with a time scale T 1 and 
an oscillatory term — cos^ ip sin^ 0) cos(|tZ)|f)/2 that de¬ 

cays with a time scale 7’2. The oscillating feature has been 
observed experimentally^. When the nuclear spin quantiza¬ 
tion axis Cz is parallel to the initial state polarization direction 
e^, i.e., 0 = 7r/2 and <p -0, the oscillatory feature disappears. 

At room temperature, when the magnetic field is along the 
z axis, the optical transition is spin conserving. The cyclic 
transition between the m - 0 ground state |g) = |0g) and 
excited state \e) = lO,,) does not contribute to nuclear spin 
dissipation since cjg - (Og - 0. When B deviates from the 
z axis, its transverse component By = + ByCy mixes 

the m - 0 sublevels and the m = +1 sublevels, so that 
(Og = -{2yg/Dgs)Rj ■ Ag and ojg = -{2yg/Dgs)RT ■ Ag, where 
Dgs (Des) is the zero-field splitting in the NV ground (excited) 
state and yg - 28.025 MHz/mT is the gyromagnetic ratio 
of the NV electron. This can be understood as a hyperfine 
enhancement of the nuclear spin g-factor— (see the next sec¬ 
tion for more detailed discussion). As a result, the nuclear 
spin dissipation rates r^(,,r± are proportional to |BtP, as ob¬ 
served experimental! y2ii^. Eor the '^Cb nucleus studied by 
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FIG. 1. Sketch of the seven energy levels of the NV center under 
optical pumping. 


Dutt et al^, the HFI tensor has been obtained by ab initio 
calculations^^ as 


KC^Ch) ^ 


-8 

0 

-0.7 




0 

-8.99 

0 


MHz, 

(8a) 

-0.7 

0 

-8.00 




-3.78 

0.19 

-1.47] 



0.19 

-5.83 

0.22 


MHz, 

(8b) 

-1.47 

0.22 

-4.12 




From this FIFI tensor, we estimate |rjg - rje| ~ 0.3 MHz when 
|Wg| = 1 MHz. Under optical pumping rate R - y\ (the exper¬ 
imental condition^), the two-level fluctuator model [Eq. (I4al l1 
gives a nuclear spin dephasing rate l/7’2 ~ F^ » (250 /ts)“* 
[« 1 /tnv ~ (4 ns)“', so the NV center is a good Markovian 
bath], which is two orders of magnitudes smaller than the ex¬ 
perimentally observed value ~ (1 Equivalently, to be 

consistent with the experiment^, the difference \(Og- a)e\ must 
be assumed to be 10 times larger—. This large discrepancy 
suggests that the leakage from m - 0 subspace tom - ±1 sub¬ 
space may plays an important role in determining the nuclear 
spin dissipation. In the next section, we shall demonstrate 
that the small leakage from |0) subspace to the | + 1) subspace 
could introduce additional contributions that may dominates 
the nuclear spin dissipation. 


III. SEVEN-LEVEL ELUCTUATOR MODEL: 

NUMERICAL AND ANALYTICAL RESULTS 

A. Model 

Now we consider general optical pumping of the NV cen¬ 
ter at room temperature, incorporating the finite intersystem 
crossing between m - 0 and m = +1 subspaces. In this case, 
the electron-nuclear coupled system still obeys Eq. O. The 
only difference is that now the Liouville superoperator X.e for 
the optically pumped NV center includes seven energy levels: 
the ground triplet \m)\g) = \mg) (i.e., |0g), | + 1^)), the ex¬ 
cited triplet \m)\e) = \me) (i.e., |0e), | + U)), and the metastable 
singlet \S}, where |g) (|e)) denote the ground (excited) orbital. 
The unitary part of X.e is a seven-level NV Hamiltonian (How¬ 
ever, we haven’t consider the effect of ionization of the NV 


center which would happen when the laser intensity is very 
strongi^) 


where A is the optical detuning between the zero-phonon 
line and the laser frequency, Hgs = ■ Sg and 

Hes - + TcB ■ Sg describe, respectively, the ground 

state triplet with zero-field splitting Dgs = 2.87 GHz and the 
excited state triplet with zero-field splitting Des = 1.41 GHz. 
The dissipative part of Xg includes various dissipation pro¬ 
cesses between the seven levels of the NV center as sketched 
in Eig. [T] the spontaneous emission from the excited orbital 
\e) to the ground orbital |g) with rate ji - 1/(12 ns)2S, the 
non-radiative decay from | + Ig) to the metastable singlet \S} 
with rate « yi followed by the non-radiative decay from 
\S) to |0g) with rate y^ - 1/(143 ns)2S, the leakage from |0g) 
to I + Ig) with equal rates ys 2 7ji, and the orbital dephasing 
of the excited state L - &e,e with rate ~ 10^ MHz2i21. 

As discussed in the previous section, there are two pro¬ 
cesses contributing to nuclear spin dissipation. The one in¬ 
volving the electron spin flip is strongly suppressed away from 
the NV center ground state and excited state level anticross¬ 
ing. Thus, in our analytical derivation below, we consider 
the other process that does not change the electron spin state, 
i.e., we drop the off-diagonal electron spin flip terms in F 
and only keep the diagonal part. The magnetic field com¬ 
ponent Bt = B^Cx -t- ByCy perpendicular to the N-V axis (z 
axis) slightly shifts the electron levels and mixes the electron 
states from \mg} and Iwg) to \mg) and Iwg) {m = 0,+1). Eor 
yglBxl Dgs, Des, the level shift can be safely neglected, but 
the state mixing has a nontrivial influence on the diagonal part 
of F, i.e., we need to keep the terms diagonal in the mixed ba¬ 
sis Ithg) and |mg): 


'y ] (^mg,mg{fhg|Sg|/lig) ■ Ag -H i7,^{irZg|Sg|irZg) ■ Ag). 


Up to the first order of the small quantities |ygBT|/£>gs and 
|ygBT|/£>es, we obtain (hereafter \mg/e) stands for |mg/g)): 

(OgISglOg) - -^Bt, 

<0g|Sg|0g) ^ -^Bt, 

^es 

<+lg|Sg| + Ig) ss +ej -H 

Uga 

(±U|Sg| + Ig) ~ +ej + Bj. 

^es 

The terms proportional to Bx lead to hyperfine enhancement 
of the nuclear spin g-factor-, e.g., the term 0 j( 0 g|Sg| 0 g) ■ Ag 
in F can be written as d- 0 g, 0 jyvBx ■ [-2ygAg/(yAfDgs)], where 
[■ • ■ ] is the correction to the nuclear gyromagnetic ratio by the 
HFI conditioned on the electron being in |0g). Since yg|Bx| ^ 
Dgs,Des, the hyperfine enhancement terms in (+lg|Sg| + Ig) 
and (+lg|Sg| + Ig) can be safely dropped, so that 

F ~ iSg^j^g E B e,z^e tr 0^,0^ Itg + CTo^^Ogltg, 
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where = -{lyelD^^)^^ ■ A.g, = -{lyelD^^)^^ ■ A.e and 

bg/e = Cj • ^gje- 

Now the second term of Eq. ([T]) describes the nuclear spin 
precession conditioned on the electron state: the precession 
frequency is yN^ + ( 7 a:B + be) when the electron state is 

I + 1 j) (I + le)), or yivB ± (Ta^B + a^) when the electron state 
is |0j) (|0e)), or y^vB when the electron is in the metastable 
singlet |5). The optically induced hopping of the electron 
between different states randomizes the precession frequency 
and leads to nuclear spin dissipation^^. Below we derive an¬ 
alytically a closed equation of motion of the nuclear spin to 
describe these effects. 


B. Lindblad master equation for nuclear spin 


The time scale tnv = 1 / min{yi 2 ,^) for the seven-level 
NV center to reach its steady state is determined by the time 
scale of the slowest process: the intersystem crossing from 
m = 0 subspace to m = +1 subspaces if the optical pumping 
is strong, or the optical pumping rate R = from the 

ground orbital to the excited orbital if the optical pumping is 
weak. When tnv is much shorter than the time scale Ti, ^2 of 
the nuclear spin dissipation, we can follow exactly the same 
procedures as used in the previous section to derive the Lind¬ 
blad master equation for the nuclear spin density matrix and 
the Bloch equation for the average nuclear spin angular mo¬ 
mentum. The former (latter) has exactly the same form as Eq. 
(I2]l [Eq. O] and describes the nuclear spin dissipation in the 
tilted cartesian coordinate (ex, Cy, Cz) with Cz = w/|t[)| [see 
Eq. Q] and (b - yArB + Po^Hg + Po^ag, where Pq^ and Pq, 
steady state populations of the NV center on |0^) and lO,,). The 
detail expression of steady populations is given in appendix A. 

Now we discuss the analytical expressions for the nuclear 
spin pure dephasing rate T,^ and relaxation rate r±. The for¬ 
mer comes from the fluctuation of Fz, while the latter comes 
from the fluctuation of F± = Fx ± iFy- Since F is a lin¬ 
ear combination of Sg^^,Se,z, and (To^,o,, the fluctuation 

of Fz and F± involve various cross-correlations among these 
four operators. Eortunately, due to the large orbital dephas¬ 
ing at room temperature, the optical pumping rate from the 
ground orbital to the excited orbital is nearly independent of 
the spin state and the coherence between electron states can 
be neglected. This allows us to neglect the cross correlation 
between (5^,^, and [(TOj.Oj, (see Appendix B for de¬ 

tails). So r,^ and r± are the sum of the contributions T^*', T^'^ 
from the fluctuation of Sg^z’Se,z associated with the m - +1 
subspace and the contributions T^' from the fluctuation 
of o’Oj.Oj, associated with the m = 0 subspace. Unless 
explicitly specified, hereafter we consider a typical situation 
\(b\ <S R, y,i. 


The contribution from m - ±l subspace is 


r(i) — 


-(1) 


2P_ 


7si 

P-y 

Jsl 


(be.i 


+ b 


yi +ysi+R 


g,z- ^ 

yi -H y^i -H 
R 


bg,zK^z— 


I 

b,,x) -(b,,yb,,y)^ 


( 10 ) 


where P, is the steady-state population of the electron state |i). 
Eormally tJjO and T^^ are independent of the magnetic field, 
but actually the components bg^Zi bg,j. s bg^ex + bgjCy, etc. 
are defined in the tilted coordinate ex, Cy, ez [see Eqs. 0], 
which in turn depends on the magnetic field. Importantly, eJ,*' 
and rO' (Jo not vanish even in zero magnetic field. This pro¬ 
vides a possible solution to the puzzling observation of nu¬ 
clear spin dephasing in zero magnetic fielc^, which has been 
speculated to be due to the orbital fluctuation of the NV center 
in the excited statai^. 

Equations (|9l) and (ITOl i exhibit four features. First, E^^' and 
do not vanish when = hg and = a^, as opposed 
to the two-level fluctuator model [Eq. 0]. This is because 
in the two-level fluctuator model, the electron only hops be¬ 
tween Ig) (with nuclear spin precession frequency y^vB -H a^) 
and \e) (with nuclear spin precession frequency yA(B -H a^): 
when ag = a^, the electron hopping does not randomize the 
nuclear spin precession, so there is no nuclear spin dissipation. 
By contrast, in the seven-level fluctuator model, the electron 
can hop between seven energy levels, each of which corre¬ 
sponds to a different nuclear spin precession frequency (see 
the discussion at the end of the previous subsection). There¬ 
fore, even if the hyperfine of the excited state is the same as 
that of the ground state, the dissipation process also exists. 
This conclusion is different from the expectationi^ii^ that the 
decoherence comes from the hyperfine difference between the 
ground state and excited state. The nuclear spin dissipation 
vanishes only when all these precession frequencies are equal, 
i.e., when = bg = be = 0. Second,r^^^ and E^^^ are 

proportional to the electron population cx y ^2 in 

the I + le) level, which vanishes when the leakage rate ys 2 
from m - 0 subspace to m = +1 subspace vanishes. Third, 
under weak pumping R « yi,yii, we have E^^^E^'^ cx IjR 
increasing with decreasing pumping strength, until the pump¬ 
ing is too weak for the Markovian assumption tnv ^ Ti,T 2 , 
based on which our analytical formula are derived, to remain 
valid. Upon further decrease of the pumping strength, the NV 
center becomes a non-Markovian bath and the nuclear spin 
dissipation rates would show a maximum and then decrease 
(see the next subsection for more discussions). Finally, under 
saturated optical pumping, E^*^ and E^*^ are saturated instead 
of being suppressed: 


p(i) 


-d) 


2x ^ 
2T 

2x5. 

4T 


^g,Z + De,Z 


(lla) 

(llb) 


where t = Ijy^x -H1 / ys 2 +1 /y^ ~ 1 Iys 2 is the average duration 
of one electron hopping cycle, fi s 2/yji is the uncertainty 
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of the dwell time in the | + 1)(| - 1)) level, and the prefactor 
2 accounts for the contribution from the m = +1 and m - -I 
subspaces. 

Equations (fTTli can be understood as follows. First, under 
strong pumping, the hopping time between the ground orbital 
and the excited orbital is negligibly small, so the m = +1 (or 
m - subspace effectively becomes a single energy level 
with nuclear spin precession frequency + (b^ + be)/2(or 
ywB - (bg + be)/2). Second, the duration t x l/yj 2 of one 
hopping cycle is ultimately limited by the slowest process: the 
intersystem crossing from m -Qiom - ±\ subspace. There¬ 
fore, Eqs. (fTTT i correspond to an effective two-level fluctuator 
model [cf. Eqs. Q]: one state is the m = H-l(orm = -l) sub¬ 
space with nuclear spin precession frequency 'yArB-i-(b^-i-bc)/2 
(or y^rB - (bg -H be)/2) and the other state is the subspace ex¬ 
cluding 1 + 1) subspace, which produce a nuclear spin preces¬ 
sion frequency yvrB. Although strong optical pumping sup¬ 
presses the randomization of the nuclear spin precession due 
to spin-conserving electron hopping between the ground or¬ 
bital and excited orbital inside the m - +1 (or m = -1) sub¬ 
space, there is extra contribution due to the random electron 
hopping between the m - ±\ subspace and the subspace ex¬ 
cluding 1+1) subspace. 

The contributions from m - 0 subspace involve and 
Ae in a quadratic form, so E^^^E^^ cx |BtP increases signif¬ 
icantly with the magnetic field components perpendicular to 
the N-V axis. Due to the finite leakage from m = 0 into 
m - ±\ subspace, the analytical expressions for E® and E^^ 
are very tedious (see Appendix B), so here we discuss the 
limits of weak pumping and strong pumping. Under weak 
pumping, oc IjR decrease with increasing pumping 

strength (this behavior does not persists down Xo R ^ 1 /Ti or 
1 /7’2, where our Markovian assumption does not hold). Under 
strong optical pumping, they are saturated: 


r(0) ^ Ii 

iP 2f 

p(0) _ Jo_ 
± 4f 


'Ag +ae 

V 2 



(12a) 

(12b) 




FIG. 2. Comparison of analytical (solid lines) and exact numerical 
results (dashed lines) for nuclear spin l/r 2 [(a) and (c)] and l/Ti 
[(b) and (d)] in a magnetic field B, = 5 mT along the N-V axis as 
functions of the optical pumping rate R. Relevant parameters are 
Ag/i; = Ag/cC^Ci,)lr] (j] - 1, 10, 100), andy .,2 = y.5i/25 in (a) and (b); 
Ag/e = Ag/f('^Ci,)/10 and yi 2 = 0,'ysi/W00,y,i/25 in (c) and (d). 


we have E 


■d) 


= E*'^ = 0 and 


T-(0) 2 

Uu - - (ag,z-ae,z) , 


-(0) 


2R + yi 
1 Po,Pog 
22R + y\ 


(a„,i - a,,i)A 


(13a) 

(13b) 


where Pq, = 1 - POg - RI(2R + yi). This recover the room- 
temperature two-level fluctuator model [Eqs. (HI and (|5|]. 
This can be easily understood: since the population is trapped 
in the m - 0 subspace, the fluctuation of the nuclear spin 
precession frequency could only come from the difference be¬ 
tween ag and Sg. 

Below we discuss two situations: (i) the magnetic field is 
along the N-V axis (z axis); (ii) the magnetic field is perpen¬ 
dicular to the N-V axis (z axis). 


where 

V2 [2 i ~y.,«ys, 1 

ffl — —~ y -trr ~ — 

JsiT \yfi ysijs y, 

is the uncertainty of the time for the electron dwelling at the 
m - 0 subspace. Similar to the contributions from the m - 
+1 subspace, under strong optical pumping, the contributions 
from the m - 0 subspace correspond to an effective two-level 
fluctuator model: one is the m - 0 subspace with nuclear 
spin precession frequency y^B + (ag + He) 12, the other state is 
the subspace excluding the m - 0 subspace, which produce a 
nuclear spin precession frequency y^B. 

When the leakage from m - 0 subspace to m = +1 subspace 
is neglected (i.e., ys 2 — 0), the steady-state populations in 
the m = +1 subspace vanish, corresponding to perfect optical 
initialization of the NV center into the state |0g). In this case. 


C. Magnetic field along N-V axis (z axis) 

When the magnetic field is along the N-V symmetric axis 
(z axis), we have ag = a^. = 0, so the average precession 
frequency th = y^yB is along the -z axis, and the tilted axis 
(ex, Cy, ez) can be chosen as (e^, -Cy, -e^)- Since ag = a^. = 0, 
only m - ±1 subspace contribute to nuclear spin dissipation: 
E^ = E^^^ and r± = E^^^ [see Eqs. (0 and ([l0|]. 

To begin with, we demonstrate the validity of our analytical 
formula Eqs. 0 and (ITOl by comparing them with the exact 
numerical results from directly solving the electron-nuclear 
coupled equations of motion [Eq. 0]. we estimate the typical 
nuclear spin dissipation time ~ Ti,T 2 ^ tnv for In this 
case, the NV center is a highly non-Markovian bath beyond 
the description of our analytical formula. To see how our an¬ 
alytical fomula becomes progressively applicable when going 
from the non-Markovian regime to the Markovian regime, we 
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Magnetic field(mT) Magnetic field(mT) 


FIG. 3. Magnetic field dependence of the nuclear spin studied by 
Dreau et alJ3.: (a) Ti and (b) T 2 times from our analytical for¬ 
mula (solid lines) and numerical simulations including (black dashed 
lines) or excluding (orange dashed lines) the electron spin-flip terms 
in F. The experimental results (empty squares) is also shown for 
comparison. The two arrows indicate the ground state and excited 
state anticrossings. 


manually scale down Ag and A^, by a factor 77 = 1 , 10 , and 100 
to decrease the nuclear spin dissipation. Figuresl^a) and|2tb) 
show three features; (i) Both the exact results (dashed lines) 
and our analytical results (solid lines) tend to saturate at large 
R, even for very strong HFI {rj - 1), where the NV center 
is highly non-Markovian. (ii) With increasing 77 and/or opti¬ 
cal pumping rate R, the nuclear spin dissipation rates l/ri ,2 
decrease and/or the electron dissipation rate 1 /tnv increases, 
thus our analytical results begin to agree with the exact nu¬ 
merical results, (iii) For successively small R, the analyti¬ 
cal dissipation rates l/Fi ^ (solid lines) tend to diverge, while 
the numerical results (dashed lines) exhibit a maximum value 
~ 1 /tnv- This is because at sufficiently small R, the time scale 
of the NV dissipation tnv ~ 1 //? is longer than the nuclear 
spin dissipation and the NV center becomes a non-Markovian 
bath. In this case, the electron-induced nuclear spin dissipa¬ 
tion rates 1 /Fi and 1 / 7’2 are upper limited by the electron dis¬ 
sipation rate ~ 1 /tnv- 

In Figs. |2lc) and |3d), both the exact numerical results 
(dashed lines) and our analytical formula (solid lines) show 
that the nuclear spin dissipation rates 1 /Fi ,2 increase rapidly 
with increasing leakage rate 7^2 from 777 = 0 to 777 = +! sub¬ 
spaces, due to the rapid increase of the population [see 

Eqs. (|9]l and ([T0l) 1. When jsi = 0, the population = 0, 

so our analytical formula gives vanishing nuclear spin dissipa¬ 
tion rates, while the exact numerical results give a extremely 
small dissipation rates. This residue dissipation comes from 
the process involving the electron spin flip, which have been 
neglected in our theory since it is strongly suppressed by the 
large electron-nuclear energy mismatch away from the ground 
state and excited state anticrossing. Nevertheness, for ex¬ 
tremely small ys 2 (= Tii/lOOO), it is responsible for the small 
difference between the analytical results (blue solid line) and 
the exact numerical results (blue dashed line) at large optical 
pumping rate in Fig. |2ld). 

Next we study the magnetic field dependence of Ti and 
T 2 and compare them with the experimental measurements^^. 
For the nuclear spin dephasing time ^2 in Fig. lUa), away 
from the ground state and excited anticrossing of the NV cen¬ 
ter [indicated by arrows in Fig. [Bja)], our analytical formula 
agree well with the exact numerical results, whether or not the 


electron spin flip terms in F is included. This indicates that the 
contribution involving the electron spin flip is negligibly small 
compared with the contribution not involving the electron spin 
flip. 

In deriving Eq. ([TOl i for the nuclear spin relaxation rates, 
we have neglected the small nuclear spin level splitting | 7 a;B|. 
When this effect is included, the analytical expressions for F± 
are given in Eq. (IA2I) . which shows a Lorentzian dependence 
on the magnetic field F± oc l/(|Bp -H bg) with a characteristic 
width 


^ ^ I_ ^ _2iL 

^ ^{2R + yi)^+2{R + yi)ys\+y]^\7N\ 

Under saturated pumping, as is usually used for optical read¬ 
out, this width ~ ti/Ita'I ~ 500 mT. By contrast, the contri¬ 
bution involving the electron spin flip also has a Lorentzian 
dependence on the magnetic field, but with a much smaller 
characteristic width ~ y\lye ~ 1 mT. The magnetic depen¬ 
dence of the relaxation time of nucleus has been measured 
by Dreau et ah^. They found that the anisotropic components 
^g,zx - Ag,xz Ag^^y - Agj^ of the ground state HEI signif¬ 
icantly contribute to the nuclear spin relaxation. The compo- 
ment ;. has been measured to be 0.25 MHz, while the other 
components are not clear. Here we assume A^, = Ag with an 
isotropic diagonal components = ‘^g,y,y - ^g.z.z - 0-25 
MHz and a small anisotropic component Ag^^^x = ^g,x,z =1-5 
kHz and A^ ^j, = Agj,^ = 0. Eigure[3[b) shows that the exact nu¬ 
merical results obtained by directly solving Eq. ([1]) agree well 
with the experimentally measured Ti time^i. As discussed 
previously, the exact numerical results contain two contribu¬ 
tions; the one not involving the electron spin flip (which is 
treated by our analytical formula) and the one involving the 
electron spin flip (which is not treated by our analytical for¬ 
mula). Eigure[3b) shows that our analytical formula provides 
an accurate description to the former contribution, although 
in the present case the latter contribution dominates because 
of the much larger isotropic HEI ~ 0.25 MHz compared with 
the anisotropic HEI ~ 1 kHz. Einally, for the nuclear spin 
at lattice O as reported in the supplement of Ref. it has 
a much shorter relaxation time ~ 40 ms at 200 mT. Such 
short relaxation time is obviously dominated by the mecha¬ 
nism of Eq.dTOll ,from which, we can estimate the anisotropic 
HEI component of this nuclear spin to be ~ 20 kHz. 


D. Magnetic field perpendicular to N-V axis 

Without loosing generality, we consider the magnetic field 
B = ByCy along the y axis of the conventional coordinate. In 
this case, the precession frequencies = —{2yeByjDgs)^y- Ag 
and ae — -{2yeBylD^^)ty ■ Ag are proportional to the magn- 
tic field. The nuclear spin precession frequency cb = y^vB - 1 - 
Pa^ag + Po,ae deviates from the z axis. In this case, both 
E^^^F^^ [see Eqs. (|9]l and ( [T0i )1 from the m - +1 subspace 
and F^°\F^^ from the m - 0 subspace are nonzero. Eor F^^^ 
and F^*\ the quantities bg^z, bg,j_, etc. are defined in the tilted 
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Optical pumping rate R (MHz) 



External magnetic field By (mT) Evolution time (ps) 



FIG. 4. Comparison of analytical (solid lines) and exact numeri¬ 
cal results (dashed lines) for nuclear spin IIT 2 [(a), (c)] and l/Ti 
[(b), (d)] in a magnetic field By = 10 mT along the y axis as 
functions of the optical pumping rate R. Relevant parameters are 
Ag/i; = AgieC^Cb)lr} (rj = 1 , 10 , 100 ), and 7,2 = ~/sil25 in (a) and (b); 
Ag/g = Ag/gC^Ci,)/lO and 7,2 = 0, 7.51 /1000, 7^1 /25 in (c) and (d). 


FIG. 5. (a) Nuclear spin 1/Ti (solid line) and I/T 2 (dashed lines) 
from numerically solving Eq. O compared with the experimentally 
measured decay for the initial state being an eigenstate of (squares) 
and /j (triangles), (b)-(d) show the dissipative evolution under (b) 
By = 0.1 mT, (c) By - 2 mT, and (d) By = 10 mT for the initial state 
being an eigenstate of Tj (solid line), ly (dashed line), and (dotted 
line). The parameters are R = 6.4 MHz and 7^2 = 7 . 51 /30. 


coordinate ex,eF, Cz = w/|w| that differs from the conven¬ 
tional coordinate (e^, Cy, Cj). 

First, we compare our analytical formula for the nuclear 
spin l/Ti and l/ 7 ’ 2 tothe exact numerical results from directly 
solving the electron-nuclear coupled equations of motion [Eq. 
O]. To see how our analytical fomula becomes progressively 
applicable when going from the non-Markovian regime to the 
Markovian regime, we start from the strongly coupled nu¬ 
clear spin and downscale its HFI tensors Ag('^Ci) and 
AeC^Cb) [see Eqs. ([^] by a factor rj - 1, 10, and 100 to de¬ 
crease the nuclear spin dissipation. The nuclear spin IIT 2 and 
1 IT 1 shown in Fig. |4]show very similar behaviors to the case 
when the magnetic field is along the N-V axis [cf. Fig. (|2]l], 
including the saturation at large optical pumping rate R and 
the improved agreement between the analytical results and the 
numerical results with increasing 77 and/or R. In particular. 
Figs. Uc) and|4jd) show that IIT \^2 increase rapidly with the 
leakage rate 7 j 2 , indicating that in addition to the contribu¬ 
tions and from m - ±l subspace, the contributions 
r® and from the m - 0 subspace also increase with 7 j 2 . 
For 7^2 = 0, the nuclear spin dissipation becomes very slow. 
In this case, the contribution from the processes involving the 
electron spin flip (not included in our analytical treatment) is 
no longer negligible. This leads to the discrepancy between 
the analytical results (black solid lines) and the numerical re¬ 
sults (black dotted lines) in Fig. |3d). 

Finally we set the scale factor tj - I and compare our 
theoretical results with the experimental measurements^. In 
this case, the strong HFI makes the NV center a highly non- 
Markovian bath, so our analytical theory only provides a qual¬ 
itative description for the nuclear spin dissipation. Since 
Ag(^^Cb) and Ad^^Ch) [see Eqs. ([8])] are approximately 
isotropic, and are almost along the y axis, while and 
be are approximately along the Cj. axis. For relatively large 


By, the magnetic field term ybiByCy and the HFI contribution 
a^, Se oc ByCy dominates the average nuclear spin precession 
frequency d>, so the nuclear spin quantization axis Cz oc ru 
is almost along the y axis. Since a^ and a^ (b^ and b^) are 
nearly parallel (perpendicular) to Cz, the m = 0 (;« = + 1 ) 
subspace mainly contribute to the nuclear spin pure dephas¬ 
ing (relaxation), so that T^^ » F® increase quadratically with 
the magnetic field, while F± » F^^* is nearly independent of 
the magnetic field. In other words, we expect that the nu¬ 
clear spin 1/7’2 to increase appreciably with By and the nu¬ 
clear spin 1/Ti to be nearly independent of By, as confirmed 
in Fig. 13 a). According to the Bloch equation Eq. (|7]), since 
the experimentally used initial states are eigenstates of 4 and 
Ix, their decay time is largely determined by 72- Indeed, for 
By » 1 mT, Fig. 13 a) shows reasonable agreement between 
the numerically calculated I/T 2 and the experimentally mea¬ 
sured decay time of different initial states. Note that the two¬ 
fold degenerate I/T 2 correspond to identical decay of (lx) and 
(4) [see Fig. |3d)]. By contrast, for By —> 0, the average 
nuclear spin precession frequency tj is dominated by a small 
term (5g_j.)bg - 1 - {Se,z}be along the z axis (neglected in our ana¬ 
lytical treatment). In this case, the fluctuation of b^ and bg of 
the m - ±l subspace mainly contribute to nuclear spin pure 
dephasing, while the fluctuation of and of the m = 0 
subspace mainly contribute to nuclear spin relaxation. Corre¬ 
spondingly, in Fig. |3a), the nuclear spin relaxation l/Ti oc 
vanishes at By = 0 , while the nuclear spin l/ 7’2 is two-fold de¬ 
generate, corresponding to near identical decay of (lx) and (/, ) 
[see Fig. 13^)]. Due to the switch of the nuclear spin quanti¬ 
zation axis at intermediate magnetic field By ~ 1 mT, the as¬ 
sociation of the solid line with 1 /T 1 and the dashed lines with 
1 /7’2 in Fig. |3a) near the crossover region is meaningless. 
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IV. CONCLUSION 

We have presented a numerical and analytical study for the 
nuclear spin dephasing and relaxation induced by an opti¬ 
cally illuminated NV center at room temperature. When the 
NV center undergoes a single cyclic transitions, our analyti¬ 
cal results provide a physically transparent interpretation that 
substantiates the previous resultsi^ and demonstrate the pos¬ 
sibility to control the nuclear spin dissipation by tuning the 
magnetic fieldi^. For general optical illumination of the NV 
center incorporating finite inter-system crossing, our numeri¬ 
cal results agree with the experimental measurements^^. Our 
analytical results suggests that the random hopping between 
the m = 0 (or m = +1) triplet states and the corresponding 
remained subspace of the NV center could significantly con¬ 
tribute to nuclear spin dissipation. This means that increasing 
the spin polarization degree of NV center would effectively 
suppress the optical induced dissipation process. This contri¬ 
bution referred here is not suppressed under saturated optical 
pumping and provides a possible solution to the puzzling ob¬ 
servation of nuclear spin dephasing in zero magnetic field^. 


Appendix A: Lindblad master equation of nuclear spin 


{F)e = POg^g + Po,^e for the seven-level model. Under this 
decomposition, Eq. ([T]l becomes 


p^£eP-i[io-lp]-i[f-lpl (Al) 

where d> = y^vB -H {F)^ is the total magnetic field that defines 
the nuclear spin quantization axis. Consequently, the nuclear 
spin dephasing and relaxation should be defined in the carte¬ 
sian frame (e^, Cy, Cz), where Cz = tI)/|o>|. 

Third, we decompose F ■ I into the sum of the longitu- 
dional part Fziz and the transverse part (F+/- -H F_/+)/2, 
where 0± = Ox ± iOy- Then treating F ■ I by the adiabatic 
approximation^^ up to the second order gives Eq. (|2]i for the 
nuclear spin density matrix p(t) = Trep(f), where 




= Re r 
Jo 


Tr,Fz{eP^'FzP)dt = Re<Fz; Fz)o 


is the nuclear spin pure dephasing rate due to the fluctuation 
of Fz at zero frequency, and 

1 1 
= - Re = - Rt{F^- 

^ Jo ^ 

is the nuclear spin-flip rates due to the fluctuation of F^f. at the 
nuclear spin precession frequency |d;|, and 


Here we derive a closed equation of motion for the nuclear 
spin from the coupled equation of motion Eq. ([TJ, where £e is 
the Liouville superoperator of the two-level or seven-level NV 
model. First, we calculate the steady state density matrix P of 
the NV center from £eP - 0. For the two-level model, the 
steady state populations on \e) and |g) are Pg - RI{2R + ji) 
and Fg = 1 - P^, where R = is the 

optical transition rate from |g) to \e) and 6^P{x) - {jln)l{x^ + 
y^) is the broadened (5-function. For the seven-level model 
at room temperature, due to the large orbital dephasing rate 
y^p ~ 10^ MHz, the spin-conserving optical transition rates 
from the ground orbital |g) to the excited orbital \e) are all 
equal to R OpRly^ for different spin states. The steady-state 
population on |0g) is 

p _ R + y\+ 2yg2 _ 

2R + ri+2r,2(2^i^ + f)' 


The populations on other NV levels are 
R 


Pog 

P±h 


-Fo., 


F -I- 7i -H 2ys2 - 

R 7^2 

- - P±lg ^ —POg, 

R + 7i+7si 7si 


and Pg ^ (2.7sil7g)P±ig. When the leakage from the m - 0 
subspace to the m = +1 subspaces are neglected by setting 
7 s 2 = 0, we have F±i^ = F±i^ = Rj = 0 and Rq, = 1 - Pog = 
F/(2F -H 7 i), which recovers the two-level fluctuator model. 

Second, we decompose the HFI into the mean-field part 
(F)e ■ I and the fluctuation part (F - (F)^) ■ I s F ■ I, where 
(F)g s TrFR is the average Knight field from the NV cen¬ 
ter, e.g., {F)^ = Pg^g + F(,tt)e for the two-level model and 


J r*+oo 

l'rede^'^‘~“^^'bPdt - -Txed{£e - i(F)~^bP 

0 

is the steady-state correlation at frequency a> between the fluc¬ 
tuation a = d - Tr dP and the fluctuation b = b - Tr^R. 

For the seven-level NV model, we have 


Fz - bg^zS g,z + beyS e,z + ^g.zd"Og.Oj + fle.zd" 

F± ~ bg ±Sg^^ -{- bg ±Se,z ^g,±^0g,Qg^ 

where O = O - Tr OP is the fluctuation part of electron op¬ 
erator O, ag^± = Ggy ± iagj and bg^± = bgy ± ibgj, etc. We 
can verify that the group Se,z, §g^z and the group (To^.o,, d’Og.Og 
have vanishing cross-correlation, so F± = F^*^ -i- F^' and 
Fp = F^*^ + F^°* can be written as the sum of the contributions 
from the m = +1 subspaces; 

F^') = Re(|V|2<5,,,;5,,,)o + \be,z\HS epj e,z)o) 

+ Re bgybey ^{Sg^z^ Re,z}o {Re,z^ Rg,z}o^ 5 

~ e,z'-,S e^z')±\oi\) 

1 . . . . 

g,z-’ ^ e,z}±\(o\ ^g,T^e,±{S S 

and the contributions from the m - 0 subspaces: 

F® = Re((3[^_2('^o„oJ<3-o„o,)o + ag,z{o‘0j,0j; o-0j,0j)o) 

+ Re%zae,z({o-o„o,;d-0j,0g)o + (d'Og.Og', d-Og,Og)o), 

~ 2R^{\^e,A.?'{d'Qg,0g',&0g,Qg)±\ib\ + |a^,x|^<(T0j.,0jd-0j.,0j>±|5)|) 

+ 2 R6('2g,=FQe,±(o'0„0j O'0j,0j)±|a)| + agyOe^Ti^ Oj.Oj O’0.,0,)±|a)|)- 






10 


Using the analytical expressions for the above correlation 
functions in Appendix B, we obtain as Eq. (| 9 ]l and 


rL‘’ = —/(|w|)|b,,j^(i + 


r.|2 


^ + yi + 7s\ 1^1 
Rjs 


R 


(A 2 ) 


+ ^/(IruDIbj,^ 
7^1 


R 


r>|2 


P-U r....... aP + yi+Tsi R + yi+jsi 


R 


Rjsi 


jrn-wru u ■.2R + 2ji+y,i 

■I-/(l<2^l)(bg,± ■ bc,±)--- 

ysi ^ 

+ -/(l<<^l)(bg X be)z— - - - 

yj'i ^ ^ 


1 - 


Rys 


where 


/(«) 




R^y,i + [(2R + 7i)^ + 2.{R + 7i)7s1 + 


Appendix B: Steady-state correlation functions 

Here we use the equation of motion method to evaluate the 
eight correlation functions. For example, the correlation func¬ 
tion <d-o„o.;(To.,o,)o can be written as -TreO-o,,o,A = -Ao„o,, 

I 


where X = -Ce^d-Q^fl^P obeys Tr^A = 0 and Xij = {i\X\j). The 
large orbital dephasing rate 7^ ~ 10 ^ MHz allows us to ne¬ 
glect the off-diagonal coherence of the electron and only keep 
the diagonal populations E, s The equations of mo¬ 

tion of Xij is obtained by taking the (/, j) matrix element of 
X.eX = d-Q^fl^P. We find that the equations of motion of the di¬ 
agonal (off-diagonal) elements of X involves the off-diagonal 
(diagonal) elements. By eliminating the off-diagonal elements 
in favor of the diagonal elements, we obtain 


-(71 -H 2yj2 + R)Xo^fi, + R2^0g,0g - Eo,(l - PoJ, 


ysXs,s + (71 + R)Xo^fi^ - RXo^flg - -Eo^Eoj, 

-ysXs,s +7ii("^-k,-k +A+4,+4) = -PsPq^, 
“(71 + 7j1 + R)X-\^,^i, -h= -P-yPog, 
“(71 + 7j1 +^)A+i,,+i, -h= -P+i^Pog, 

(71 + -1, + 7 j 2 Ao„o, = -P-igPo,^ 

(71 + + 7 i 2 Ao„o, = -P+igPo,, 


where Ajj is the energy difference between the electron 
state |/) and \j} in the rotating frame of the pumping laser. 
Solving the above equations gives the correlation function 
{o‘o„o,;o-o„o,)o = as 


, PoPo, 1 PoX^-Po,-PoP ( R + js 1 1\ Po.(E-i^ -H P+i^) 1 PoPs 1 

(cro„0,;cro„o,)o - -1-;- -1-+ -;- — -;-■ 

l+ri 2 R + ji 1 - 1 -77 \2f;-l-7i7., 7^1/ 1 -1-77 2 R + yi I + J] y,i 

where 77 = 27^2/7^1 + 2 (/? -1-275)7^2/(75(2/? -1-71)) is a dimensionless constant much smaller than unity since 7,2 <si 7si, 75. Using 
the same method, the other correlation functions are obtained as: 


(O'05,05(0O'05,05(0))0 - 


-I- 


PogPo, 1 

fi - 2752//? 1 n) 

, , POgil - POg - POg) 1 

fl , 7.1 1 

1 /? + 71 + 2752 

1 -H 77 ' 

[ 2/?+ 71 Rj 

' 1+77 ' 

[R ' 75 ( 2 /? + 7 i)J 

7*1 


P OgiP-i, + P+iP 1 R +yi + 2 ys 2 POgPs l/?-i-7i-i- 2 ys 2 
1 + j] R 2 R + yi 1 + T] R 751 


<o'o„o.(f)o- 05 , 05 ( 0 )>o 


POgPOg 1 ^ Pogjl - Pq, - Pog) I 1 ^ R 1 \ Po,(P-i, -Hf+i^) _ Pq^Ps 1 

1 -I-77 2 /?-I- 7 i ’’’ 1 - 1 -77 \75i ’’’2//-H7175/’'’ ( 1 -i- 77 )( 2 /?h- 7 i) 1 -h 77 7si ’ 


(o‘o.,o,(0o'o„o.(0))o 


Pog(l-Pog)l+2yg2lR JOgPogi 

{ 7*2 1 7*2 '1 

, , P0g(P-\g + P+\g) 1 

fl l-275i//?\ 

1+77 2 /? + 7 i 1 + 77 ' 

\Rys\ ( 2 /?+ 71 ) 75 ] 

1+77 1 

\R 2/?+ 71 / 


PoXl-Po 


1 + 7 ] 


■ Pog) (R + yi ^ R + yi 1 


Rysi 2R + yiygj 1 + rj 


POgPs R + yi + 2ys2 


Rysi 


Similarly, the correlation functions at finite frequency are ob¬ 
tained as 


{R e,z ^ R e,z}(i} — 


2 (R + ito) 


Rysi + i(y\ + 7ii + 2/?)w - 07 


:P-U, 


A , 2 (/? -H 7 i -I- 751 -H ico) 

{ g.z, g,z)^ /?75i +7(2/?+ 71+751)07-^2 -h’ 

. . _ 2 R 

{Re,z,Rg,z}o, +i( 2 R +yi -t-7ji)w- 

2 (/? + 7 i) ^ ^ 

Rysi + i( 2 R + yi + ysi)co - to^ 


{R g,z '^ R e,z}oj — 


When the leakage from the m - 0 subspace to the 777 = +1 
subspace are neglected by setting 752 = 0 (and hence rj - 0), 
only the populations Pq^ and Pq^ are nonzero, so only the first 
term in the above expressions survives: 


<o'o„o.;d-o„o.)o = <o' 05 , 05 (?)o' 05 , 05 ( 0 ))o = -<d' 05 , 05 (f)d-o.,o,( 0 ))o 


- -<o-o,,o,( 0 o‘ 05 , 0 j( 0 ))o 


P^P% 

2 R + yP 
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which give Eqs. (fTJt of the main text. For saturated pumping, and hence Eqs. (fTSl i of the main text, 
we have 

(o' Oj.O*; O' o„o,)o = (o'0„0,; O' o„o,)o 
= (o'0„0,;0'0j,0j)o 


- (o'Og.Oj; o'o„o,)o 

-1 

8f 
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